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In a previous paper, we showed that the problem of tunnehng in quantum wires was integrable 
in the isotropic case ga ~ 2. In the present work, we continue the exploration of the general 
phase diagram by looking for other integrable cases. Specifically, we discuss in details the manifold 
5p + So- = 2, where the associated "double sine- Gordon" model is integrable. Transport properties 
are exactly computed. Surprisingly, the IR fixed points, while having complete reflection of charge 
and spin currents, do not correspond to two separate leads. Their main characteristic is that they 
are approached along irrelevant operators of dimension 1 + j- and 1 + -j- , corresponding to transfer 
of one electron charge but no spin, or one spin 1/2 but no charge. 



I. INTRODUCTION. 



While much progress has been made recently in the field of quantum impurity problems, it is still an open challenge 
to fully understand the phase diagram for a single (charge) impurity in a one dimensional quantum wire, where the 
electrons are described by a Luttinger liquid with charge and spin degrees of freedom. 

The problem was studied using the renormalisation group and perturbation theory by Kane and Fisher and 
Furusaki and Nagasoa ||]. A complete picture of the flows generated by the RG was obtained in terms of the "g" 
factors, gp and describing the interaction of bulk charge and spin degrees of freedom. The existence of tentalizing 
new fixed points, which would be partly transmitting and partly refiecting, was in particular conjectured. While steps 
were taken to identify these fixed points using conformal field theory |^ , the determination of their exact nature - in 
particular the computation of the universal associated conductances - remains open. 

In this paper, we are addressing this impurity problem from the point of view of integrability. A simpler problem, 
the case of an impurity in a spinless electron gas, was mapped on the boundary sine-gordon model Q, which is 
a completely integrable quantum field theory In our case, the problem of electrons with spin and impurity, it 
maps after bosonization and the same folding maneuvers as on a "double sine-Gordon model" , with two coupling 
constants Pi and (32 (related simply to gp and g^). Unfortunately, this model does not seem to be generically integrable. 
There are however special integrable manifolds. In the first part of this work |^ the case g^ — 2 was solved by using 
a slightly different theory than the double sine-Gordon model, with however the same boundary interaction. Here a 
more direct route is adopted and wc investigate the integrability of the "double sine-Gordon models" . The simplest 
one is given by gp + ga = 2, and can be considered as a particular case of a new integrable model recently discovered 
by Fateev . 

This integrable manifold is, a priori, not the most interesting, since it lies far away from the region where the new 
IR fixed points are conjectured to arise. However, a closer study along the lines of Q gives rise to very intriguing 
results. By a variety of arguments, including the computation of the non equilibrium tunneling current at vanishing 
temperature, we found that the IR fixed point does not correspond to two separate leads ("open" fixed point [^|), 
as suggested by a (maybe too naive) analysis of the classical action . It is not completely clear to us what this IR 
fixed point is. Although it corresponds to complete reflection (vanishing transmitted current) for spin and charge, 
the operators along which it is approached seem to correspond to collective excitations that would be impossible for 
separate leads: the leading operators describe actually either transfer of one electron charge but no spin, or one spin 
1/2 but no charge! 

The paper is organized as follows. Section 2 is a brief summary of the field theoretic description of the problem and 
of the folding of Q. Section 3 studies the "double sine-Gordon model" in full generality. The integrable manifold is 
studied using the technique of non local conserved currents of |^, and S and R matrices are obtained. The IR fixed 
point is (partially) identified and shown to differ from the expected one. In section 4, the out of equilibrium current 
with an external potential is computed exactly. Section 5 contains our conclusions. Many additional technical details 
are given in appendices. In appendix A, we explain why the model, although quantum integrable, is presumably not 
classically integrable. This is in sharp contrast with the case of the ordinary sine-Gordon model, and suggests that 
semi classical arguments as in Q are maybe not reliable. Appendix B discusses the problem of quantum boundary 
integrability. In appendix C, a perturbative computation of the S matrix around the free fermion point is presented, 
which confirms our non perturbative identification based on non local currents. Appendix D is a first order Keldysh 
computation that fixes the relation between the bare coupling in our action and the parameter Tb of the R matrix. 
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Appendix E sketches the computation of the whole two point current correlator using form- factors, which allows a 
more complete identification of the operator content of the IR fixed point. 



II. THE MODEL. 

In order to fix notations and the problem, let us repeat some of the definitions in In the one dimensional 
Luttinger liquid, bosonisation of the fermionic operators is accomplished via 

n odd 

with /i =1, j. The fields and 0^ have commutation relations 

[^^ix),ef,,ix')]^iS^^^,0{x-x') (2) 

from which we can devise two different representations of the Luttinger liquid (our conventions follow |^). Here, we 
will work in the so called 9-representation. Changing basis to the charge and spin degrees of freedom 

Gp^et + ei, e^-Gt-ex, (3) 

we obtain the action 

S^Sp + S^ (4) 
= / dxdy^iid^Qpf + {dyQpf] + -L[(5,e,)2 + [dyO^f]. 

where y is the imaginary time, the spin and charge velocities have been normalised such that ~ Vp ~ I. In this 
convention, the g-factors have value ga — gp — 2 for a non-interacting system. 

The electric and magnetic conductances of this system follow directly from Kubo's formula 

Gp/a- = "^5p/<T- (5) 

As discussed in [||, this formula is actually not correct for the charge conductance, because of the way the system is 
coupled to the reservoirs. One rather has 

Gp = 2^, (6) 

for any gp. In this paper, we compute the conductance as in ^ to make comparison with the results of this paper 
more straightforward. The correct physical results follow by simple rescaling. 

We are interested in the physical properties of this model in the presence of a charge impurity at the origin a; = 0. 
Due to this impurity, the hamiltonian gets an additional piece 

SH^ J dxV{x){ij\,p^ + V-jVi). (7) 

Here V{x) is a potential which has essentially zero measure outside a: = and under the bosonisation rules sketched 
at the beginning, this leads to the change of action 

SS = X J dy cos ^/nQp cos \^Qcr, (8) 

where A ~ V{2kp). In general, all terms consistent with the symmetries of the action must be considered at the 
boundary since they will be generated by renormalization. This means that the most general form of the perturbation 
is 

dS^ Jdy J2 ^^e^-^("-«''+""«') (9) 
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where the A's are real couplings. The renormalisation group equations read at first order 



In this paper, for reasons given later, we will restrict to the manifold 

5^+5p = 2, (11) 

so that only the perturbation with Hp — ria = ^ is relevant. In the following we denote An by A. 

In [Q, the equivalent problem for spinless electrons was exactly solved. The solution required a folding to transform 
the impurity problem into a boundary problem, and then used recent results on boundary integrable quantum field 
theories |^ together with the massless scattering approach. The same folding can easily be acomplished in the problem 
with spin. First, it is convenient to rescale the fields, writing the action as 

S — Sp + Scr + \ J dy COS y/ng^Qp COS y/Wg^Q a- (12) 

s^ = I dxdy^iid^e^f + {dye^f]. (13) 

We then introduce odd and even fields 

= i= [Qp^Lix, y) + e^,fl(-x, y)] (14) 

% = ^[Qp^L{x,y) ^ QpA-^,y)]. 

With this, the interaction at a; = involves only the even fields which have Neumann boundary conditions. The odd 
fields, having Dirichlet boundary conditions, completely decouple and do not interact. In the previous form, the even 
field is left moving and we can "fold" using 

Ot,,L = %{x,y), X < 

6pM^%{-x,y), x<0, (15) 
and now expressing everything in terms of the fields 6 which is defined on the negative axis, we get the action 

S = So + B 

- f dxj dy\ J2 [(^-^m)' + idy(^,)'] + ^ J dy cos ^J^Op{0) cos sj^O^iQ). (16) 

The next step in Q was to use the integrability of the corresponding boundary quantum field theory - the boundary 
sine-Gordon model (BSG). In the present case with spin, the problem, involving two fields, is more complicated, and 
very few results are available. We must first investigate in more details this "double sine-Gordon model" (DSG) and 
its boundary counterpart (BDSG). The following section, together with a series of appendices, are devoted to that 
study. The reader interested in physical properties only can go directly to the end of next section. 



III. THE DOUBLE SINE-GORDON MODEL 

A. Non local currents and integrable manifold 

Although the problem at hand involves an interaction only at the boundary, it is crucial - like for the ordinary 
sine-Gordon model |]l^ - to understand first the model with a bulk interaction. We start thus with the general action 
for the double sine-Gordon model, with two fields (pi, (f>2 

S ^ Sq + Sfjuik 

= 5*0 + A / dxdy cos I3i4>i cos [32(t>2- (17) 
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A more general form of this action was first investigated in [Q, where an integrable manifold was identified, of which 
our model is a subset. The following arguments, in particular the use of affine quantum group symmetries, are ours, 

but were strongly influenced by the results of . The double cosine term (the perturbation) has dimension x = ^^^f^ ■ 
The free propagators are {4>i{z, z)(j)j{z', z')) — —4^% In \z — z'p, where we set 2: = y — ix. 

Unfortunately, we do not expect the model to be integrable for any choice of the coupling constants /9i , /?2 ■ A useful 
strategy to identify quantum integrable varieties is to look for non-local conserved currents, following Bernard and 
Leclair A current J{z) will be conserved if it is local with respect to the perturbation, and if the residue of the 
simple pole in the short distance expansion of J and the perturbation is a total derivative. By trial and error, the 
simplest possible candidates have the form 

Jf(z)=5^,; e^^t^^ (18) 

where ^Pi{0) are the chiral and antichiral components of the field, and in ([l8|), i ^ j ■ One has for instance (exponentials 
are implicitely normal ordered) 



z — w 



47r(z — w) 



-ij^V2{z) ^ill3iipi(w)+l32'P2(w)] ^ ^i[l3lipiiw)+l32V2{w)-^ip2{z)] 

so the simple pole has residue 
This is a total derivative iff — ^ = /Jt — Ir- or 

fil + ai = 4^. (19) 

In that case, the perturbation has physical dimension d = 1 Q. Setting 

djf =dHf- djf =mf, (20) 

the currents give rise to conserved charges 

Qf cx / dzjf + / dmf 



f (X y dzjf + y dzHf. (21) 

By counting arguments, following [pr] , one expects that the conservation extends to all orders in perturbation theory 
^ . In addition, we of course have two conserved topological charges ^ 

T, = ^ r d,<p. (22) 

J-00 

The invariance of the action constrains Ti and T2 to have the same parity. 

Following the construction of ||], one checks that the non local conserved charges generate the quantum afhne 
algebras slq.{2), with 

f). 



^In the case of the ordinary sine-Gordon model, this dimension corresponds to the free fermion point. The manifold ( p^ does 
not, however, correspond to free fermions except when /df = 2tt. 
^The fact that the dimension of the perturbation is a rational number weakens the argument of 
^Note that we divide by vr, not 2tt as in the sine-Gordon model. 
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Moreover, one also checks that charges of different type i ^ j commute. For instance, one has 



so the simple pole has residue 



I3i{z - w) 



dipi 



^ 4:Tt\ dip2 , dipi 



0i) Pi ' 



P2iz - W) 



'ft 



which is a total derivative on the manifold (|l9|). 

Again, while the commutation relations are established to first order in perturbation theory, simple scaling argu- 
ments show that they hold to any order. We conclude that, to any order in perturbation, the theory exhibits the 
symmetry slq^ (2) (g) slq^ ( 2) 

Under the condition (tt9[), it is lengthy but straightforward to show the existence of non trivial local conserved 
quantities (this has first been shown in H). Besides the stress energy tensor, the next non trivial one has dimension 
4, and follows from 8X4 = 9i?4, where 



^ alii - SaDidipiY 



4at af(5Vi) + (1 -> 2) + 6aiaiidipidip2) 



(24) 



anda^^^. ... . . . „ 

The existence of this non trivial conserved quantity is expected to be enough, following standard arguments |12[ to 
guarantee quantum integrability. We can make an easy guess for the particle excitations and the factorized S matrix. 
Commutation with the two quantum afhne algebras independently dictates a tensor product structure. Moreover, 
since Ti , T2 must always have the same parity, excitations should be massive particles carrying a pair of quantum 
numbers Ti,T2: ie they are simultaneously solitons of the field (jfi and (1)2 . One then obtains as the simplest choice 



S = -5s 



'5- 



where 8/3 is the S-matrix of solitons for the ordinary sine-Gordon model at coupling /3, 

. 2 



(25) 



(26) 



and we used the results of This scattering matrix was first found in [Q using a different argument. It is important 
to stress that we do not have independently solitons of the field 4>i and of the field (/)2, but composite objects. Also, 
the coupling in the S-matrices is renormalized non trivially. This is largely the consequence of the dtp term appearing 
in the conserved currents, a feature which is not there in the standard sine-Gordon model, but is necessary here 
because the interaction is the product of two cosines: although the S-matrix is factorized, the theory is thus highly 
interacting. It is interesting to discuss the action and the factorized S matrix ( p5| ) from the point of view of a fermion 
theory: see appendix C for more details. 

Some special points can be identified on the integrable manifold, which can be solved independently. For instance 
the case Pi = (32 = V^n can be recast, by forming linear combinations of the fields, into a problem of two decoupled 
sine-Gordon theories at their free fermion point, in agreement with (p5|). The case /3i = is also trivially solvable, and 
equivalent to one free boson and one sine-Gordon model at the free fermion point. This is a very different description 

from the one provided by (|2^), which involves in particular a sine-Gordon S matrix at the point (32 — "3^, ie the 
N — 2 supersymmetric point. Nevertheless, we will check that, when physical quantities are computed, they have 
identical expressions in both descriptions. 
From the relation 

A 

27r 



Ti,e 



(27) 



we find the soliton operators (up to polynomials in derivatives of 1 







exp 



exp 



(28) 



which have topological charge Ti = ±1. The local current jf' is a raising (lowering) operator for particles of type i: 
it annihilates solitons (antisolitons) and transforms antisolitons (solitons) into solitons (antisolitons). jf' conserves 
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B. Boundary integrability and new IR fixed points 



The main step in jj] was to use the integrabihty of the boundary sine-Gordon model. Here, one can easily prove 
that the model ( [l7| ) defined on the half line x G [—00, 0] with a boundary interaction at a: = 



B = X 



dy cos Y0i(2/) cos ^(l>2{.y)-, 



(29) 



is integrable for any choice A, A. This model we call the double boundary sine-Gordon model (DBSG). In the factorized 
scattering picture, the boundary interaction is described by a reflection matrix solution of the boundary Yang Baxter 
equation. The factorized structure of the bulk S-matrix leads to the immediate conjecture for this reflection matrix 



R 



R 



P2 



(30) 



where R is the reflection matrix of the boundary sine-Gordon model given explicitely in 

Of course the original Kane-Fisher impurity problem is massless in the bulk. As in the spinless case, we describe 
this problem starting from the massive theory with a bulk coupling ([T^) by using the quasiparticles basis, and letting 
A ^ 0. In that limit, we get right and left moving massless solitons with dispersion relation e — ±p = me^ , with 
the rapidity. Movers of opposite types interact via a constant S-matrix, while movers of the same type still interact 
with (|2^). In addition, R movers bounce back as L movers at the boundary. We recall in particular the physical 
amplitudes 



(%)= 



1 



^ 1 



(31) 



Here, Tb is an energy scale related with the coupling A by c>c A^, where we used (|l|). 

Already at this point, the identification of these amplitudes allows a rough analysis of the infrared fixed point. 
If we expand the modulus square of the R-matrices in terms of the inverse coupling we find that l-R^P expands in 



powers of Tg'^^^^' oc A ^^/ft . On the other hand, by the Fermi Golden rule, these are proportional to the square 
of the matrix elements of the operators perturbing the IR fixed point. We thus see that these operators must have 
dimension 



277 



1. 



(32) 



The same result is obtained when analyzing physical properties (see below) . 

This conclusion is, a priori, very surprising. Indeed, the standard analysis of the IR limit based on the classical 
action (|9|) would suggest that, at large coupling, the fields (pi are located at the bottom of the boundary potentials, 
leading to Dirichlet boundary conditions with 



2n 

i{x = 0) = —^ni, rii + rtj = mod 2. 

Pi 



(33) 



This fixed point we refer to in the following as the standard IR fixed point. With (|33|), it is easy to classify the allowed 
operators describing the approach to the standard IR fixed point. One simply computes the partition function of the 
model with these boundary conditions at a; = 0,2: = L, and matches with the formula of conformal invariant field 
theories, Z — Trw^°~'^^'^^. Here, one finds 



Z = 



1 



ri{i 



E 



1 +— s 



(34) 



where the primed sum indicates that rii,n2 have the same parity, w = e ^/^'^, The dimensions of the allowed 
operators, up to integers, are thus of the form 



1 are 



In particular, the lowest dimension operators have dimension || -f- Remarkably, the ones with x = ^ 



not there! 
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There are several explanations for this discrepancy. The first, of course, is the possibility that our S and R 
matrices have not been correctly identified. The previous arguments about affine quantum group symmetries are 
however usually considered to be quite strong. They can be made even stronger at some special points with N — 2 
supersymmetry (see below). Moreover, physical results for all limiting cases where the solutions is known by other 
means will be reproduced correctly (see also below). Finally, the S and R matrices also agree with perturbation theory 
around the free fermion point j3f — 2tt (see appendix C). If we do believe these S and R matrices, another possible 
explanation is that the perturbation around the IR fixed point has a vanishing radius of convergence, so the expansion 
of (^l|) (or other physical quantity) in powers of Tb is the result of a resummation, and cannot be used to naively 
identify the dimension of the perturbing operator. Such phenomenon has been observed in other integrable models 

We suspect however that this discrepancy arises, more simply, because the IR fixed point obtained by naive analysis 
of the action is not the correct one. To justify this, we recall first that integrable flows are usually found to be integrable 
considered either as perturbations of the UV or of the IR fixed point, and with the same set of conserved quantities. 
For instance, in the boundary sine-Gordon model perturbed by cos/9</>, the local conserved quantities like commute 
not only with the perturbation, but also with its "dual" cos ^4>- As a result, the fiow near the IR fixed point, which 
is approached along this dual field plus local integrals of motion T2n, is integrable too. In the present case, we checked 
that our T4 quantity does commute, besides cos/3i(/)i cos/32</>2, with the currents , but not with cos ^^i, nor with 

cos ^'/'i cos ^02- As a result, it seems impossible to flow to the disconnected fixed point, since its approach would not 
be integrable. On the other hand, the IR fixed point described above was found to be approached along operators of 
the same dimension as jf^ . Very presumably, it is approached along these operators indeed (plus the T2„), providing 
an integrable fiow. 

As another piece of justification, it is useful to consider a generalization of (^7|) given by the action 
S=\j dxdy [(9^</)i)' + {d^^hf + [d^^f] dxdy {e"^ cos(/3i0i + M2) + e'"-^ cos{Mi - fS^h)] (35) 

involving thus one more bosonic field. Non local conserved currents generalizing (fiSh can be built 



provided the condition 



J±(z) = (iP.d^, ± f3d^) e^'^'"^' (36) 



Pf + Pi -P^ = 47r, (37) 
holds. By the same arguments, the model is integrable. An integrable boundary action follows 
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(38) 



It is possible to compute the algebra generated by the conserved charges and it is found to be the same as in the 
previous model. Thus, the S and R matrices for this model are essentially the same as the ones for the previous, 
simpler model with /3 = 0, with however the the new condition (^). The IR fixed point deduced from the classical 
action is the "standard" one like for /3 = since the field cannot be compactified. 

This generalization is useful because at the value P\ = Att, the model ( ^ ) coincides with the bosonized form of 
the N = 2 supersymmetric sine-Gordon model, as can easily be inferred from the fact that the Ji currents have 
dimension 3/2: = G±. The bulk perturbation is unitary and preserves supersymmetry [|T^ . The same is true 
for the boundary perturbation. This supersymmetry is a particular case of the quantum affine symmetry generated 
by the currents . The preservation of supersymmetry is a more firmly based property than the preservation of 
quantum affine symmetries however, because the generators are local. We now use the conservation of supersymmetry 
to justify the identification of the IR fixed point provided by the S and R matrices. Let us first recall what happens 
for the supersymmetric point of the ordinary boundary sine- Gordon model first [|5). Then, /J^ = 1|ZL. The perturbing 
operator cos/3<^ in the UV has boundary dimension x ~ 2/3 = ^ + ^ corresponding to the superpartner of the 
chiral primary field X, and preserves supersymmetry. The IR fixed point is, in that case, correctly reproduced by 
the standard analysis of the classical action. The Dirichlet boundary conditions analogous to (^) give rise to the 
partition function 

^=^E-'^=XoT, (39) 
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which is the Neveu Schwartz character of the identity (that one is in the NS sector follows immediately from the fact 
that boundary conditions are the same on both ends of the cylinder). The IR fixed point is approached along the 
operator cos^t^ of dimension x — 3/2, ie along the operators G^. This is very natural from a Landau Ginzburg 

point of view: near the minima of the potential, the effective potential is W = with no left-over chiral primary 
field in the spectrum beside the identity. The approach has thus to be along a descendent of the identity, with 
as the natural candidates. It is not totally obvious a priori that a perturbation along G^ preserves supersymmetry; 
while G"*" anticommutes with G^, the anticommutator of G'^ and G~ is non zero. However, the G^ are fermionic 
operators; to have a bosonic IR action, they must thus be coupled to some fermionic boundary degrees of freedom, 
and these allow the restoration of supersymmetry [p^ . 

We can now get back to the model (^) when Pf = An, ie with unbroken N = 2 supersymmetry. It is easy to check 
that the standard IR fixed point (p3[) cannot be approached while preserving supersymmetry. On the other hand, 



the leading operators describing the approach to our new IR fixed point have dimensions d = 3/2 and d — 1 + 

The first one corresponds to = G^, and, like for the ordinary sine-Gordon model, it does preserve supersymmetry 
after addition of fermionic boundary degrees of freedom. The second one corresponds to . Because the two affine 
quantum algebras in the problem commute, it too does preserve supersymmetry. Hence, our fixed point is perfectly 
compatible with unbroken supersymmetry, while the standard one is not. 

Away from the supersymmetric points, a similar argument can be made based on the conservation of affine quantum 
symmetries instead |^. It leads right away to the claim that the IR fixed point should be approached along the non local 
conserved currents: J the ordinary sine-Gordon case, as is well known, and along jf, of dimensions 

di = 1 + in the double sine-Gordon case, in agreement with the conjectured S and R matrices . 

Based on the analysis of physical properties (see below and appendix E), we conjecture that the operator content 
of the IR fixed point is a subset of the spectrum encoded in 

^2 _2 

1 , ^ -27r-^ . , -27r 



2. .E- 

^ ^ ill 



w ^ (40) 



Note that the partition function ( pO| ) would correspond to a model with Dirichlet boundary conditions of the type 
( ^ ) but with no coupling between rii and n2. A possible conjecture for the ratio of boundary entropies in the UV 
and IR is 



9IB. _( Pi 0: 



guv \l3i + 2tt Pi + 2t: 



1/2 



(41) 

while the same ratio for the standard fixed point would read ( -^^^ = 2 (^^\ . If this result is correct, 

^ ) standard J 

our new fixed point is less stable than the open (separate leads) one: therefore, it might well be that it is reached 
only with the UV action (||) on the integrable manifold gp + go- = 2. 

The question of course remains: why is the identification of the IR fixed point misleading in our case, while it was 
not for the boundary sine-Gordon model? Wc have no satisfactory answer to that question at the present time. But 
there is an interesting indication: in sharp distinction with the ordinary sine-Gordon model, the double sine-Gordon 
model (|^ seems to not be integrable classically, as discussed in the appendix A. Integrability in the present case is 
thus a truly quantum property. Since integrability is also responsible (via the analysis of non local conserved currents 
and the resulting S and R matrices) of the unexpected nature of the IR fixed point, it is maybe not so surprising that 
this IR fixed point cannot be obtained from a classical analysis of the action, and thus differs from the standard one. 

C. Summary 

We can now get back to the physical problem. A comparison of (^6|) with the results of this section leads to the 
identification of the integrable manifold in terms of the g-factors: gp + ga = 2. The postulated exact S'-matrix has 
the form of a tensor product of spin and charge degrees of freedom 



*This is a little tricky to do in practice, because of non locality problems for general values of (3i. It is however possible for 
special values of f3i 
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S = ^Sp(E)S^, (42) 

with Sp_cr the 5- matrices of a sine-Gordon model with parameters ^ = ^^rf- Similarly the boundary reflection matrix 
is given by 

R = ~Rp(E)Ra. (43) 

The elementary excitations of the folded theory are quasiparticles carrying a charge number Qp = ±1 (in units where 
e = 1) and a spin number = ±1 (in units where a physical electron has = ±h). Here, we have set 



1 1 



dx dxOp 



Qa = I dx dj, (44) 



Having a complete description of the kinematics of the model, we can now compute physical properties. We will 
restrict to the T — Q case here. Extensions to finite T are in principle possible, but quite technical. 



IV. DC CONDUCTANCE AT T = 

We will exactly follow the same line of argument as in The quasiparticles are considered as a gas interacting 
through factorized scattering. Standard thermodynamics arguments can be applied, in the presence of electric and 
magnetic chemical potentials, to compute the densities and the filling fractions of these quasiparticles. These are 
expressed via the thermodynamic Bethe ansatz. Having these densities, a rate equation can then be written to 
compute the contribution to the conductances of the current backscattered by the impurity. 



A. General computation. 

We only consider the case of an external voltage, so the quasiparticles have chemical potential e^^/^^, T the 
temperature. At T = 0, a Fermi sea is formed, consisting of positively charged quasi-particles, with either spin up 
or down. These quasiparticles carry charge labels (1,1) and (1,-1). Due to the factorized form of the S matrix, 
the scattering of these quasiparticles is factorized into the scattering of the first and second labels. The first labels 
scatter diagonally with the element (■S'p)]]]^. The second labels in general scatter non diagonally. They would scatter 
diagonally if was of the form \ / integer , but then bound states (carrying a non vanishing electrical charge) would 
have to be taken into account. It is in fact simpler to analyze the more general case corresponding to model (|3^), 
taking gp — pi, g^ = P2, Pi integers. The final results can then be analytically continued to the region of interest 
( |l9|) with /3 = (a computation right at /3 = is also possible, although more complicated. We will describe this 
elsewhere). The scattering of the magnetic quantum numbers is then non diagonal. The TEA is slightly intricate to 
write: one needs to introduce pseudoparticles to diagonalize this scattering. The final result fortunately is very much 
like in the ordinary sine-Gordon model [pT^,n8| , and is conveniently encoded in the TEA diagram 



-Pl + l O P2-1 

— o- - -®- - -o — o 



-Pl O O P2 



The TEA equations are 



k 



V k 

e^Pi = 0- e-Pi+i = 2" 2^ * (45) 
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In these equations, the star indicates a convolution of the functions with the kernel k{9) = 2 corfi e ' ^^'^ ^~ designates 
the negative part of e [p^ . It is a simple exercise in Fourier transform to extract from ( ^5| ) an equation satisfied by 
eo itself. Setting e = — Cq , one has 

<e) - J^^ He - e')eie')de' ^^~\e' (46) 

Here, the kernel $ has fourier transform 

l> = P(i^_p,+2+Xp,_2) (47) 

where k{x) — 2 cosh x ' kernel ii'_pj^+2 arises from the left part of the diagram, and coincides with the derivative of 
the logarithm of (S'p)ii, and the kernel arises from the right part of the diagram 



sinh(pi - l)a; cosh(p2 - 2)a; 

' —r-r , ifp2 -2 = 2 cosh a; — — 

smhpix cosh(p2 ~ 



/f-p,H-2^2cosha; °^^^2u^ ^p.-2 = 2 cosh.T^^^^^^^ (48) 



leading to 



$ = 1- -oMp^+P--^)^ . (49) 
2 sinhpix cosh(p2 — 1)2: 



We recall that e describes the excitation energy of new quasiparticles which do not have a well defined magnetic 
quantum number - this is because the scattering of magnetic numbers is non diagonal. The total density of these 
quasiparticles is obeys 

Pie)- j\{e'Me-e')de' = l-e', (50) 

where we have reinstated Planck's constant. This density follows simply from the excitation energy as p = 

In the previous equations, A is the Fermi rapidity defined by e{A) — 0; it is the edge of the Fermi sea. A 
can be determined, and the foregoing equations solved, by using a standard Wiener Hopf analysis. If we write 
K{e) = (5(61) - $(6'), we have 



^ AT( ^Mf V (51) 

N[uj)N{—uj) 



with 



' ^' .gp^(^c^5p/2)^(l/2 + ^w(.9,-l)/2)' ^ ' 

a function analytic in the lower half plane when we choose A = l/21og2 + (g^r — l)/21og((7o- — l)/2 + gp/2\oggp/2. 
Here the conventions for the Fourier transforms are 

— «5(c^), ~9{^) = / d0e-«5(^). (53) 

Then following exactly the same steps as in Q we find the solution for the fourier transform of the density 

2h (1 + luj) 

The excitation energy of the particles follows 

e»e-- = _,^«M + !^NH)Ni.)^ ^^^^ 



The condition e{A) = is equivalent to 
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lim u;e(itj)e-'"^ = 0, (56) 



which in turns gives the expUcit value of the Fermi rapidity, A 

In the presence of the impurity, the backscattered current is computed by using a Landauer Biittiker approach [ p^ . 
The only difference with the computation done in Q is that the quasiparticles have a magnetic degree of freedom. 
However, we are only interested here in the properties of the electric charge, so we simply sum over all possible 
channels of reflection for the magnetic quantum number. This leads to the final expression for the current 

.A 

I{V,TB) = e depie)\iRp)X{e~9B)\^. (58) 

J —oc 

This expression was derived under the assumption gp,ga integers. We expect however that it holds for any values of 
these constants, in particular for the values on the manifold gp + ga ~ 2. Indeed, physical quantities in integrable 
models have (so far) always had a smooth behaviour in terms of the coupling constants, despite the very different 
possible Bethe ansatze that depend on their arithmetic properties. Furthermore, Fateev [^| has computed some 
thermodynamic properties for the more general model (|3^). One can check on his computations, by comparison with 
perturbative results or sigma models approximations, that the results are indeed analytical in terms of the coupling 
constants. 

Let us now analyze (p8|). In the UV limit, the matrices i?^ become unity and we find the result 

J(y,TB = 0) = ^5p, (59) 

as expected ^. The full result at finite Tb is given by a double expansion, one for each reflection matrix. The modulus 
squared of the previous reflection matrices is given by 

\iRp)Xm'= ^^ \ ,,, , \iR.)ti9)f= \ (60) 

+ l + exp(-^)6' + l + exp(-^)6' 

For large voltage, we just expand the previous expression and evaluate the current at each order by suming over 
residues, we obtain the following strong coupling expansion 



n=l ^ ' 

with 

Using a similar procedure and treating the case in which Tg < Y ^ we can also expand the reflection matrices by 
splitting the integral over the rapidity, Q in two pieces. After a little algebra, we find the weak coupling expansion 

I(Y,Tb) = N(0)PP / — - — ^ <^ V -1 ; — } > (63) 

^ ' 2h ^ ^ J 27r(l + iw)| ^^ ' 2n/gp + iuj 2 sinh(7rwgp/2) J ^ ^ 

with b = (Tee""^)^/^'' . We integrate the first sum by closing the contour in the lower half plane and the second term 
by closing in the upper half plane. In the second term, the zeroes of N{lu) cancels the poles of sinh(7ra;gp/2) and only 
the poles oj — 2ni with n integer, and u; = i contributes. We get 



^To get the correct version, it suffice to rescale V as mentionned at the beginning 
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I{V, Tb) = —gp + E '^^n f ^ j + C/iTs, 



with the coefficients 



N{-i)\2n 



(64) 



2/1 y n! T[-ngp]T[\/2 ~ n{g, - l)](2n - 1) smiirngp) 

C/i = -^Bin(|gp). (65) 

Let us write the first few orders more explicitely 

mTB) = — .g. |i - 2 -(-)^ + ^r(i/2 - ,,) r(,,/2)^r(,./2)^ [^) ^ -j ^''^ 

The spin conductance in a magnetic field follows from the previous analysis simply by the interchange of g^ and gp. 
In order to compare with perturbation theory and complete the picture, we need to establish a relation between the 
boundary scale Tb and the perturbative coupling A. This is done in appendix D and we find 

Tb = „ . (67) 
2sm(f.gp) 



B. Special Cases. 



When ga ~ gp = 1 the sum (61) can be done analytically, in that case we have that 



and we find the result 



hn = (68) 

2n + 1 

m Tb) = ^-^ arctan (69) 

This result can be found by other means since by writing the action, not in the charge and spin fields basis but rather 
in the spin up and down basis, the action decouples into two Luttinger actions which do not interact with each other. 
Moreover each is at the free fermion, or Toulouse, point and the resulting conductance is well known to be of the form 
given by (|69|). The two expressions can be precisely matched using ( |67| ) together with the results of |||^. 

Another special point is = 2 and g^ = 0. The spin field then totally decouples, while the dynamics of the 
charge field is described by the ordinary boundary sine-Gordon model at the free fermion point. One thus expects 
a conductance of the form (^) (with however an overall factor of 2, and a renormalized Tb). On the other hand, 
the S matrix which we claim describes the double sine-Gordon model is not at the free fermion point in the charge 
sector, but rather at the N — 2 supersymmetric point! Nevertheless, special cancellations then occur to reproduce 
the correct physical quantities. Indeed, using the expression (|6^, we have at that particular point 



In the limit gp — > 2, this vanishes for n odd. For n even, one has 



2n + 1 \sin7r5p/2 



(71) 



^Most unfortunately, there are two misprints in the latter reference. In eq. (6.20), the left hand side should have a 2XAi 
instead of Ai. The next equation should have a term eV/2TB instead of eV/Tb in brackets. 
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This reproduces the correct result using (|67|). 

A last special point is gp = 0,ga- = 2. There, one expects the spin field to essentially decouple as the fluctuations 
of the charge field become very large. Explicit computation using ( |6^ gives rise to 



-._n,^H- ,.„,._;,._,., (^^| I (72) 



This coincides indeed with the classical limit of the spinless problem 



V. CONCLUSION 



The main result of our study is the existence of a new kind of IR fixed point in this impurity problem. This 
fixed point has a vanishing transmitted current but it does not correspond to the open (separate leads) fixed point 
as demonstrated by the analysis of the operator content. Of course, since the transmitted current vanishes, it does 
not correspond to any of the other fixed points (eg current transmitting and spin reflecting, or partially transmitting 
and partially reflecting) considered so far [QH] . The main property of our fixed point is that it is approached along 
operators which describe the transfer of one electron charge but no spin, or one spin 1/2 but no charge. This is 
clearly seen by reformulating the non local currents in terms of the physical fields; a quick identification can be 
made using the dimensions d — 1 + — and d= 1 + — (recall that for the separate leads fixed point, the approaching 

9p 9cT 

operators have d — ^ + corresponding to the transfer of a physical electron carrying a spin 1/2, and d — 
corresponding to the transfer of a pair of electrons with opposite spin) . Understanding the physical meaning of these 
properties in more details is an open challenge. A possibility is that particles with half the electron charge and half 
its spin are actually transfered: "half electrons" have actually been encountered in the double barrier case, where 
their appearance corresponds to the transfer through an "island" . To make further progress, a reliable formula for 
the boundary entropy would be most useful - technical difficulties have to be resolved before this can be obtained, 
however. 

It is not clear to us how the existence of these new fixed points changes the conclusions of [^]J^. Our analysis, 
however, indicates that an instanton approach to the study of the problem might be misleading. In particular, the 
duality property, which is the main ingredient of the lore in the spinless case (although only partially understood 
yet), cannot always hold in the case with spin: the UV and IR expansions of our current do not match under the 
transformations g A/g. 

Acknowledgments We thank I. Affleck, P. Fendley, M. P. Fisher and N. Warner for useful discussions. This work 
was upported by the Packard Foundation, the National Young Investigator Program and the DOE. It was also partly 
supported by the NSF under Grant No. PHY94-07194. 

APPENDIX A: NON-INTEGRABILITY IN THE CLASSICAL CASE. 

In this appendix we want to show that the system described by the action 

r ^ 1 

5- / ^-(9^(/.,)' + Acos(/3i0i)cos(/32</'2) (Al) 

i—l 

does not have a spin 4 conserved quantity classically, and is presumably non integrable. 

Using the light-cone variables z = x + iy, z = x ~ iy, the classical equations of motion are given by 

dzdg(j)i = -;^sin(^i0i)cos(/3202), 

d.d,^2 = cos(/3i</)i) sin(/3202). (A2) 

The first non trivial condition of integrabilty is the existence of a spin 4 conserved current, T4 such that 

= 8,62. (A3) 
The most general form of the T4 current respecting the U{1) symmetries is given by 
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T4 = + ^2(9,02)" + ^3(a.</.l)2(a,02)' + Aiidl^^f + ^5(9202)2. (A4) 

When we take the dz derivative of T4, several of the terms can be ehminated up to total dz derivatives, so it is enough 
to look for a 02 of the form 

9-2 = Bi{dz(l)ifcos{Pi(j)i)cos{P2h)+B2{dz(p2fcos{f3i(j)i)cos{f32<l)2) + 

B3idz(f>i)idz(l>2) sin(/3i(/.i) sin(/3202). (A5) 

Straightforward algebra gives the relations 

-/32B2 

-P2B1 + PiB^ (A6) 
B3 
2B2. 

If Pi 7^/32, then these equations have no solutions and thus the system has no non trivial conserved quantity of spin 
4 classically. Presumably, it then does not have conserved quantities of higher spin either. 

The algebra in the quantum case is a bit different. The problem is to find a chiral quantity T4, of the most general 
form 

T4 = Aiidz^i)^ + A2(dzip2f + A:i{dzmf{dz^2f + Aiidl^if + A^{dl^2f. (A7) 

such that the residue of the simple pole in the short distance expansion with the perturbation is a total derivative. 
Using the propagators < ipi{z)ipj{z') >— ln(z — z'), one finds that the quantity T4 as given in the text (|2j) solves 
the problem. 

That some 1 + 1 quantum field theories can be non integrable classically, but are quantum mechanically integrable 
at finite values of Planck's constant {Pi here) was observed before: see for instance [Q. Usually however, quantum 
integrability becomes "classical" once the theory is fermionized. While a similar phenomenon occurs here for the 
supersymmetric case pl| , we are not aware of a classically integrable, partially fermionized, version of our model for 
pj + Pi = 47r. 



- AiKPi = - 


'PiBi, 


~A2Ap2 = - 


As^ = -P,B2 + 


P2B3, 


^ 2 


A4 2 - 






AP.P. 
^ 2 


= ^3, 


2 - 



APPENDIX B: BOUNDARY QUANTUM INTEGRABILITY 

We show that the boundary perturbation in the action ( [l6| ) preserves quantum integrability. It is sufficient to prove 
that the boundary perturbation of its conformal limit preserves the bulk integrable manifolds. The proof follows 
Ghoshal and A.Z. Zamolodchikov argument in The unperturbed conformal theory is constituted by the two 
independent bosons (pi, (j32 in the half plane, with Neumann boundary conditions 



d_ 

dx 



such that the propagators turn out to be 



{(j)^{Xl,yl)(j)j{x2,y2)) 



J_, {Vl - y2f + {Xl - X2f 



[yi - y2f + {xi + X2f 

4^ ^ k2 



(Bl) 



(B2) 



where k is a short distance cutoff. The half- line bosons 4>i{x,y) can thus be considered as the folding of two full line 
bosons ^i{x,y), 



Mx, y) = [*'(^' y) + y)] . 



(B3) 



The full fine left and right components ^i{x,y) — ^i^L{x,y) + ^i^ii{x,y) allow the further decomposition of the 
half-line bosons in a left and a right component 



<l),{x,y) = ^tix,y) + '^t{-x,y) 



(B4) 
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where we define the even fields 



whose correlations are 



^tix, y)^^ [^^,Lix, y) + ^^A~X, y)] (B5) 



m^uy.WA^.,y.)) = iog^l^^^±fci^. (B6) 

47r K 

Conformal invariance of the half-hne model is guaranteed by the condition that no energy-momentum flux escapes 
the system across the boundary, T{x,y) — T{—x,y), where T{x,y) = —2TTj2i{dz4>i) = ^27r (^^^^(a;, j/)) , 
T(x, y) = -27T = y)f, and z = y + ix. 

This condition imphes an infinite number of higher spin conservation laws, Ts+i{x,y) = Ts+i{—x,y), where the 
spin label s is an odd integer. Therefore a boundary perturbation preserves integrability if some of the conserved 
currents satisfy 

Ts+i{x,y) - Ts+i{-x,y)\^^^j = 
The standard argument |^] consists in evaluating the left-hand side at first order in conformal perturbation theory 

A J dy' [Ts+iix,y)-Ts+i{-x,y)] a>(0,y'), (B8) 

where $(0,?/) — cos /3i/2(f>i{y) cos /32/2<j)2{y) is the perturbing operator. Therefore the OPE of the currents with the 
perturbing operator are needed. Since the half-line bosons at the boundary satisfy 

Mo,y) = 2Mo.y) (B9) 



it turns out that the computation of (BJ ) is identical to the one in the bulk for the quantity SgTs+i and the perturbing 
operator cos/3i0i(O, y) cos/32'/>2(0, y). As a result, the bulk integrable manifold gives rise to a boundary integrable 
one. 



APPENDIX C: PERTURBATIVE COMPUTATION OF THE S MATRIX 

In this appendix we give more credibility to the S matrix proposed in the text by checking that it is correctly 
recovered in perturbation. As shown in the bosonic action of ( p^ ) can be put in a useful fermionic form by first 
introducing two auxiliary boson fields Xi 2 ~ 2v^('^i'?^i -'- /^202) whose action 

J {(^ + ^)[('9xi)' + (ax2n+(^-^)(axi)(5x2)+A(cos2V¥xi+cos2V^X2)} (CI) 
is the bosonization of the two fermions theory defined by 




The couplings are defined by 5 = ^(7^-1-7^ — 1) and gi — — -^). The integrable manifold (uM becomes 

Pl P2 Pl ^2 

gf = 2g{2g + 7r). (C3) 

The perturbative computation of the scattering amplitudes of the fermions permits their identification with the 
composite charge and spin solitons. Denoting by (p cr) with p,a ~ ± the four elementary excitations whose exact 
scattering matrix is (|25|), the identification is 

^1 = ( + + ), V^l = ( — ), ^2 = {+-), V^2 -(- + ). (C4) 
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The perturbative region 
in the exact solution (P5| 



for the two fermions theory imphes I3f ~ 27r and therefore each sine-Gordon component 



^{P1T1){P2T2)^^> ^ ^^P1P2 [t^,Pl)^<Ti<J2 (ojh) 



(C5) 



can be expanded around the respective free ferniion points Pi ^ A-k. We recall that the sine-Gordon couplings f3i 
are related to the couplings (3i of our model by Pi = 0^+2tt ■ example let us consider the following expansion 



(++)(—) 



,2tt 



27r 



1 



'Pi 



Pi 



sinh 6 



0(2) = - 



27rg 
sinh^ 6 



0(2) 



(C6) 



where we have also used the integrability condition (C3). This result has to be compared with the direct perturbative 
evaluation of the amplitude •S'^^^'^^- We will proceed with a dispersion relation method, as described in |2^ for the 
Thirring model. The two perturbative couplings g and gi are treated as independent. The integrability condition will 
be imposed at the end of the second order computation which is needed because o{g'l) still contribute to the lowest 
order in g once ( |C3|) is taken into account. The tree level amplitudes, schematically in figure 1, are determined by 
the two four fermion vertices in the lagrangian contracted on the respective external fermions. 




FIG. 1. Tree level diagram. 



Some of them are 



s 



s 



s 



ipiipi 

i)2i)_2 

i'2th 
Ipllpl 



ci'ii'2 

^1pl1p2 



1 — cosh 9 , 
-(l + *25— — ^) 



sinh^^ 



iAq — 



1 



sinh 6 
1 



sinh 9 
1 



sinh 9 



cosh 9 



-{l + i2g 



1 + cosh t 
sinh0 



(C7) 



S^'p{9)^tgi- 



sinh^ 
1 



sinhfl 



S^±i9)^-il + ^g.'-^Y 



4>l4'2 



sinh 9 ' 



We have written on the same lines the crossing symmetric pairs. The other amplitudes can be obtained by the listed 
above exploiting either the symmetry 1 <-!■ 2 or CPT symmetry. 

The real part of the second order contributions, schematically in figure 2, can be obtained by the cutting technique. 




FIG. 2. Second order diagram. 
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Actually it is not necessary to use explicitely the Cutkowski rule. One can use unitarity, as we show in the following 
example. Let's consider the amplitude S^^^^^- According to unitarity of the total S-matrix S''S = 1 and elasticity of 
2 — > 2 particles scattering the real part of the second order contribution is determined by the tree level amplitudes 



(s^'f + i){s'f''p y + s'^'f (5f^f ' + 1)* + s'^'f (s'f'f)* + s'f'f (sf^' y + 0(3) (cs) 



„ , 1 + cosh 6* , 1 
^8g2 gf +o{3). 

sinh^ 61 sinh^6i 

In total generality a second order diagram in the class of figure 2 can be parameterized as 



— — ^ [acosh^0 + &cosh6l + c] (C9) 

smh6/ 7rsinh6' 

where 9 — 6, in ~ 6, i-n when the total momentum P which flows into the diagram is P = (pi + P2), ±(pi — P2), 0, 
resp. Obviously the order g^ can be substituted by Qi and ggi depending on the kind of vertices involved in the 
diagram. The factor l/sinh0 is the usual jacobian of the total momentum conservation expressed in the rapidities 

01,02 and the mandelstam variable is s = ipi + V2Y = 2m^(l + coshS). The function r(9) — is the result of 

^ ^ ^ ' ^ ' TTsinhy 

the dispersion integral on the two-particles phase space 9(s — Am?)/ s{s — ^w?) 

The polynomial in coshfl is the effect of the contraction of external fcrmion wave functions on the numerators of 
the loop fermion propagators. Each diagram can at most be logarithmic divergent at high energy, log s ~ 0. This 
fact constraints the polynomial to be of order not greater than two. Therefore the amplitude of our example can be 
written in the form 

s^^^^{e)^^^+ (cii) 

4^-7: |r(6l) [A cosh^ 6* + B cosh 61 + C] -|-r(i7r- 61) Ucosh^ 6* -1-5 cosh 6* + (71 + r{iTT)[D cosh B + E]] + 
smht' 
■ 2 

{r[e)[Ai cosh^ 61 Pi cosh 61 Ci] + r{iTT - 9)[Ai cosh^ + Bi coshg -I- Ci] + r{iTr)[Di cosh 6* + Pi]} 



sinhfl 
+0(3). 

Some of the coefficients of the polynomials, which encode the freedom in the finite part of a regularized diagram, 
are determined by imposing i) the matching with the already computed real part (|C9|), ii) the asymptotic behaviour 
S^'-^'-{9) as 9 00 and iii) the crossing symmetry relation = S'^^^'^{ii: — 9). The result is 

Sff(^9) = ^ - -^(1 + coshfl) - -4 ^ { 2-4:7: cosh0 - p) + -^,E, + o(3). (C12) 



i'l^i^ ' sinh6' sinh^6' sinh^ 6* 7rsinh6' V sinh6' / 7rsinh6' 

The same procedure can be used for the amplitudes S^^^^ and S^^^^ which are mapped one into the other by 
crossing symmetry and whose asymptotic behaviour at high energy is S{9) — > constant. We obtain 

= - - """" + iS^'i - 

^ ' dcosh9 + e] '-^(-^+diCoshe + ei] +o{3) 



TT sinh 9 \ sinh 9 / sinh 9 \ sinh 9 

S^'f {9) = -l- ^%(1 + cosh 9) + -^4- [(1 + cosh 9f + 4] + 
'^I'^i^ ^ sinh6l^ ' sml?9^^ ' ^ sinh^ 6* 

*V f , , \ ^gl 



TT sinh 9 \ sinh 9 / tt sinh 9 \ sinh 9 



-4— TTT - dcoshe + e - -±— -—— - d, cosh0 + ei + o(3) 



Imposing now the threshold behaviour hm^i^^o S'^l'^li^^) = ~1 and lime^o('S'^^^^^ (0) + '5'^^^^^' (^)) = we obtain the 



four relations 
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E = -^, e^-A-d, Ei^-2-2di, ei^-l-di. (C14) 

We are thus left with only two undetermined constants, d and di. This freedom reflects the freedom we have in the 
definition of the coupling constants g and gi. Let us choose to define them through the asymptotic behaviour of the 
amplitudes 

hm 5^*(0) = -e-^^ (C15) 
y^+oo ^^^^ 

fixing the coefficients to c? = 0, di = 0. 

The amplitude we want to compare with the expansion of the exact results (p5| ) is therefore 

St^fiO) = ^ - + cosh^) - ^ - ^-^cosh^- ^ + o(3). (C16) 

0101 ' smhe sinh^e* sinh^ 6* 7rsmh6'smh0 TTSmhe* 



On the integrable manifold (C3) at lowest order in g we obtain 



501^^1 ^ ^ ^(2) (C17) 

'^i^i^ ^ sinh^e ^ ' ^ ^ 



which coincides with the lowest order expansion of the exact result (|Cq ). The same check can now be extended to 
the remaining amplitudes. 

APPENDIX D: KELDYSH COMPUTATION. 

In order to compute the differential conductance we need to use the Keldysh formalism since the system is driven by 
reservoirs. To do so, we will use the formulation of the model on the full line as described by the equations ([l^). The 
effect of the reservoirs can be implemented by shifting the charge fields, 0p ^ 8p + y/gp/7Ta{t)/2 with dta{t) — V. 
Under this prescription the current is evaluated by taking the functionnal derivative 

where we have used conventions in which h = e — 1. Using the Keldysh contour, C, which goes from — cx) to oo and 
then comes back (see figure 3), to expand the partition function we obtain to first order 



-oo c 

FIG. 3. Contour for the perturbative evaluation. 

i2. 



/(2)(o) = J^dt sin[|aW] (D2) 

where ^ = ± depending on the location of t, ie upper part of the contour or lower part of the contour. The functions 
P^''^ (t) are the corresponding contraction of the vertex operators time ordered on the Keldysh contour 

To this order the result can be found explicitely, it is given by 
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This is in agreement with the Bethe ansatz solution given in the bulk of the text provided we make the identification 
(putting h = 2iT and e = 1 in the TBA expressions and gp + = 2 in the previous perturbative expression) 

|A|2 = -sin(Jg,)rB. (D5) 

TT Z 



APPENDIX E: OPERATOR CONTENT OF THE IR FIXED POINT. 



To get a better knowledge of the operator content of the IR fixed point, it is useful to study correlation functions 
of various operators. This is done using form-factors ||2^. The form factors for theories with a factorized S matrix of 
the form (|2^) were analysed by Smirnov |Q for the massive theory described in the bulk of the text. These are now 
matrix elements of operators O in the space of states created by the Fadeev-Zamolodchikov operators 

(o|o(o,o)z;(0i)---z;j0„)|o) (ei) 

where now the asymptotic states are described by isotopic indices giving the two quantum numbers of the excitation: 
— (T'i,T2). Here we will concentrate on the current form factors. For these operators, the massless limit, which is 
needed to treat the boundary perturbation problem, has basically the same form has for the boundary sine-Gordon 
case discussed in Here, the theory has U{1) (g) U{1) symmetry, and we can construct currents invariant with 

respect to either of these C/(l). There is an infinity of non-zero matrix elements, or form factors but for the sake of 
the argument we need only the two solitons ones. Let us concentrate on the right moving currents. For the field (pi 
one has 

cosh 2^ (t^' — — in) smh 2^(6/' — — itt) 

and for the field 

The constant df^^ ^^'^ the function C/3i ,/32 are given in Q and will not be needed in the following. The computation 
of current-current correlator in the Luttinger model follows immediatly using the methods of | p3[ . By using the 
explicit expression written before, it is easy to check that < 0\d±4>id±(j)2\0 >= as expected. The zero temperature, 
frequency dependent, conductance would then follow from a correlator of the form < O|9+(/)i9_0i|O >, after the usual 
identification Op ^ (pi. Using the boundary state formulation, the previous form factors lead to a contribution to that 
correlator which is of the form 



where we used the notation Oib = Qi — Ob and Ob = logTs. Pg and are the sine-Gordon elements of the reflection 
matrix R are coupling g. We have P = — RZ and Q = Rt_ = R^. The frequency dependence follow from a Fourier 
transformation of the previous expression over time. We then get one single integral and the frenquency dependence 
can be recast in the reflection matrices by shifting the rapidities as was done in |2^]. The expansion around the IR 
fixed point then follows from the expansion of the reflection matrices. Doing this we obtain a multiple power series 
with the following exponents: (cj/Tb)^", [uj/Tb)^^^^^'^ and {lu/TbY'^"'^^^- Assuming the current correlator explores 
the whole operator content of the theory, this is consistent with (MG) 
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